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Abstract 



Gluon mediated exclusive hadronic decays of b quarks are studied within the 
^ ' standard model (SM) and the constrained minimally super symmetric stan- 

dard model (MSSM). For all allowed regions of the MSSM parameter space 
(A, tan/3, mo, W1/2) the penguin magnetic dipole form factor is domi- 
nant over the electric dipole and can be larger than the magnetic dipole form 
factor of the SM. However, overall the SM electric dipole decay amplitude 
F^ dominates the decay rate. The MSSM penguin contributions to the free 
quark decay rate approach the 10% level for those regions of parameter space 
close to the highest allowed values of tan (5 (~ 55) for which the gluino is light 
(rrig ~ 360 GeV) and lies within the range of the six d squark masses. In 
these regions the supersymmetric box amplitudes are negligible. The MSSM 
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phases change very little over the allowed parameter space and can lead to 

significant interference with the SM amplitudes. 
12.15.Ji, 12.60.Jv, 13.25.Hw 
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I. INTRODUCTION 



Supersymmetry (SUSY) is a highly favoured candidate theory for new physics beyond 
the standard model (SM). Of particular interest are the flavour changing neutral current 
transitions involving the quark-squark-gluino vertex and the non-removable CP-violating 
phases which arise as the renormalisation group equations (RGE) scale the physics down 
from the unification scale M\j ~ 10 16 GeV to the electroweak scale. These effects of SUSY 
have implications for rare B decays and mixings [|T]-§| and for other observables such as 
quark electric dipole moments 

Measurements of rare flavour- changing B decays provide opportunities for the discovery 
of indirect effects of SUSY ||[7| as the measured observables involve SM and SUSY processes 
occurring at the same order of perturbation theory. In contrast to the situation for B° — B° 
mixing where new physics is expected to change the magnitude of the CP-asymmetries but 
not the patterns of asymmetries predicted by the SM || , the effects of new physics in decay 
amplitudes depends on the specific processes and decay channel under consideration and, 
although small, may be detectable by comparing measurements that within the SM should 
yield the same quantity. 

The b —>■ s transition provides an opportunity to study CP violation from non-standard 
phases |J and there is significant current interest in the b — > sg penguin decay for which it 
has been argued fH| that enhancement for on-shell gluons is needed from non-SM physics 



to explain the CLEO measurement [TT|] of a large branching ratio for B — > rj + X s and the 
r)' — g — g gluon anomaly. 

For the gluon-mediated exclusive hadronic decays studied here the effects of SUSY are 
expected to be small and difficult to disentangle from the SM effects because of the large 
uncertainties associated with the SM predictions. The SM calculations involve the com- 



putation of the quark level decays b — > qq'q', calculation of the Wilson coefficients |L2[] to 
incorporate QCD corrections as the physics is renormalised down from the electroweak scale 
to the scale ra& and, finally, the calculation of hadronic matrix elements for the hadronisation 
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of the final-state quarks into particular final states, typically evaluated using the factorisa- 
tion assumption. As this last stage can introduce such large uncertainties that predicted 
SM rates for exclusive hadronic penguin decays can be in error by a factor of 2 to 3, we 
will restrict the present study to the weak scale quark level processes where any differences 
between SM and SUSY physics are more apparent. 

The most predictive of the SUSY models is the (constrained) minimally supersymmetric 
standard model (MSSM) [|13|,^] based on spontaneously broken N = 1 supergravity with flat 
Kahler metrics pfl , universal explicit soft-SUSY breaking terms at the scale Mmssm ~ Mu 
and spontaneous breaking of the SU(2) ®U(1) symmetry driven by radiative corrections. 
Such models contain two CP-violating phases 5^^ SM from the soft-SUSY breaking terms in 
addition to the usual phase Sqkm of the CKM mixing matrix. With the usual assumption 
that these SUSY phases vanish identically at the unification scale because of CP conservation 
in the SUSY breaking sector, it is claimed |fL5| ,|6|,py3l that the MSSM predictions for B° — B° 
mixing and penguin decays such as b —>■ qq'q' are very similar to those of the SM and that 
non-minimal SUSY models are needed to obtain any significant non-SM effects. An early 
study [TjJ concluded that superpenguins are small compared to ordinary penguins unless 



the gluino is very light (~ 1 GeV) and satisfies rrig <C mj. Recently Grossman and Worah 
|| have found that the gluonic penguin amplitudes for b — > sqq and b — > dqq in the effective 



SUSY model of Cohen et al [17] can be up to twice as large as the SM gluonic penguins and 
with an unknown phase. 

In this paper we revisit the question of MSSM predictions for the penguin mediated 
decays b — > qq'q' . In doing so we review in some detail the SM predictions with particular 
reference to the relative contributions of the internal u, c and t quarks to the gluon penguin 
PI ,the relative magnitudes of the various form factors and the role of the strong and weak 



phases p^j20|1 . We find, for example, that the CP violating phases for the b — ► dg and b —>■ sg 



electric form factors, which dominate the decay amplitude, have no simple relationship with 
any angle of the unitarity triangle. For the MSSM we explore the allowed regions of the 
parameter space to locate those regions which give the largest modifications to the SM 



results. In contrast to the SM, we always find the magnetic amplitude to dominate the 
electric amplitude. Also, there are large regions of the MSSM parameter space for which the 
magnetic amplitude is greater than that of the SM. The search for SUSY would be greatly 
aided if the magnetic amplitudes could be experimentally isolated. 

Conservation of the gluonic current requires the b — > qg vertex to have the structure 

r> 2 ) = f§ u q {p q )T a vM 2 )u b {p b ) (i) 

where 

+ia^[F 2 L {q 2 )P L + F 2 R {q 2 )P R ]. (2) 

Here F\ and F2 are the electric (monopole) and magnetic (dipole) form factors, q = p b — Pq 
is the gluon momentum, Pl(r) = (1 =F 7s)/2 are the chirality projection operators and 
T a (a = 1, . . . , 8) are the SU(3) C generators normalised to Tr{T a T h ) = \$ ab - 
The h — > qg vertex is 

r;(l 2 ) = -||§ v b ( Pb )T%(q 2 )v q (p q ) (3) 

where has the form @ with the form factors F± 2 R (q 2 ) replaced by Fi' 2 R {<l 2 ) where the 
relationship between the F and F form factors will be discussed later. 

To lowest order in a s the penguin amplitude for the decay process b — > q g — > q q'q' is 

M Pe„ g = [u^T^Mp^Ap^T^P?)] (4) 

where 

+ l ^f[F 2 L (q 2 )P L + F 2 R (q^P R }. (5) 



This gives the free quark decay rate 



dT p ^ _ 1 .g 2 , 2 1 T , 2 . n , 2m 2 at( 2 



dq 2 2887r 3 Mvr 2/ E b ^ /y q 2 



+ (6) 



where 



I(q 2 ) = [I + (m « ~ 4 g } " [1 - ^] 1/2 (7) 

is the phase space factor and 



m\ ml q 2 



N(q 2 ) = (q 2 p b ■ p q + 2 p b ■ q p q ■ g)(|F 1 L | 2 + \F*\ 2 ) - 3m b m q q 2 (Ff F* * + c.c.) 

-3m b p q ■ q{F^F* * + F?F% * + c.c.) + 3m q p b ■ q(F^F^ * + F^F* * + c.c.) 
+^(4 p b ■ q p q ■ q - q 2 p b ■ p q ){\F^\ 2 + |F 2 ^| 2 ) - 3m b m q (F 2 L F* * + c.c.) (8) 

with 4m 2 / < q 2 < {m b — m q ) 2 . 

Similarly, for b — > q q' q', the amplitude is 

M Pens = || ^(p 9 )^V 6 (p 6 )]« ? ,(P 9 0T^%,(p g -0] (9) 

where 7^ is obtained from (|5|) by the replacement of all F(q 2 ) form factors by F(q 2 ) form 
factors. The decay rate dT Peng /dq 2 is given by (§|)-(PD with the same replacements. 
One CP violation observable of particular interest is the partial rate asymmetry 



where 



dTs/A = l dT dT 
dq 2 2 { dq 2 dq 2 '' 1 ' 



II. THE GLUON PENGUIN IN THE STANDARD MODEL 

Although the expressions for the SM form factors are known, we believe it is useful to 
present a brief outline of their derivation, both in order to clarify their regimes of validity 
and to aid our later generalisation to include the effects of SUSY. 

For the SM the contributions to the b — > qg vertex from W and scalar exchange (Fig. 
§ give [0,^] 
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vfV) = E t[(K + 4) p L + B s ,p R ] (12) 

i=u,c,t 

where \ bq = g|^_ K* q K ib and X is the CKM matrix. For the 6 — > qg vertex V^, the CKM 
matrix elements are replaced by their complex conjugates. 

After putting the external b and q (= d, s) quarks on mass shell, V^ M has the form 

V^ M {q 2 ) = {ap btl + bq^ + C1 ,)P L 

+(d Pblt + eq ll + f'y li )P R (13) 

and the form factors for arbitrary q 2 are given by 

2{m 2 q - m 2 b )F^(q 2 ) = (a + 2b)m q + (d + 2e)m b (14) 
2(m 2 q - ml)F^{q 2 ) = (a + 2b)m b + (d + 2e)m g (15) 

and 

Fftf) = a/2, F 2 R (q 2 )=d/2. (16) 
Neglecting terms of order m^/M^ and ml/M^ then 

a + 26 = -m, a SM , d + 2e = m b (3 SM (17) 

where 



-1 rl-X 



a SM (q 2 ) = £ A? / <** / " X dy[4(x + 2y)(y- 1) 
~^ Jo Jo 



together with 



_2x t (l -x-xy- 2y 2 )]/Y z (x, y) (18) 

(3 SM (q 2 ) = V A- 9 C dx f l ~ X dy[4x 2 - 8x 
i Jo Jo 

+8y 2 -8y + Ylxy + Xi(4x 2 - 6x + Ay 2 

-8y + 10xy + 2)]/Y i (x,y) (19) 

W) =»,E ^ t dx C" dy[2x{\ - y) 
t Jo Jo 

+x i (l-x-xy)]/Y i (x,y) (20) 



and 



F 2 R (q 2 ) = m b £ A 

i 

+Xi(2x 2 - 3x + 3xy + 1)}/Yi(x, y). 



1 rl—x 

dx / dy[2x(x + i/) 
o Jo 



(21) 



In the above 



Yi(x, y) = x + Xi(l - x) + g 2 [xy + y(y - 1)]/M^ 

where = mf/M^. 

If these expressions are evaluated at g 2 = we have /3 SM = a SM and 



(22) 



F?{0) 



gi 



8M 2 , 



E K iq K ib fl( X i)> F?(0)=0 



m 



1 #(°) = =1? (o) = k*** / 2 i 



m 6 



8M 2 , 



(23) 
(24) 



where 23.24 



12(1 -x) 



-[18a; - 29a; 2 + 10a; 3 + a; 4 



f8-32x + 18aT)inxl 



—x 



4(1-2;)' 



- [2 + 3a; — 6a; 2 + x 3 + 6x In x] 



(25) 
(26) 



For small Xj, f2{%i) ~ w whereas fi(xi) ~ — |lnxj. 



For 6 — > gg, (g 2 ) max = (m& — m g ) 2 rs 20 GeV 2 and the assumption q 2 <C m 2 which would 
justify the replacement of the form factors with their values at q 2 = is invalid for Fi(q 2 ) 
for the u and c quarks. This observation has also been made in PJI8H- For these light quarks 
we can evaluate F±(q 2 ) by neglecting m 2 compared to m 2 in (|i~4[) and in the numerator 
of (H) so that 



fi(x h q 2 ) 



1 /-l-iE 

dx / <i?/[2x 2 — 4x + 6xy 
o ao 



+4y(y-l)]/Y i (x,y). 



(27) 



This integral is dominated by the logarithmic singularity near x = so we can set x = 
everywhere except in the leading term of the denominator to give 



8 



1 r l 



A(x J ,g 2 )^4/ dxl dy— v{ ] V) (28) 

Jo Jo x + Xi-q 2 y{l-y)/M^ 

10 2, 2 2 2^ + 1 , , 

lnx l + g(zi) (29) 

9 3 3zi 3 Zi yv ' K ' 



where Zi = g 2 /4m 2 and 

{\/— arctanL/^M, z < 1 
V * W^' (30) 

For g 2 > 4m 2 , g(z) becomes imaginary due to the generation of a strong phase at the uu 
and cc thresholds |T9|j20| . Our result for fi(xi, q 2 ) is equivalent to that obtained by Gerard 



and Hou [ 20| . For the u quark, Zi is large and we use the asymptotic form of (^9 



•2 



10 2 n _, g 2 



/ 1 (x U)g > i --[ln(-^)- l 4 (31) 



We will be concerned with the b — ► dg'g' and 6 — > sg'g' transitions. Although the form 
factors -F\ and — F 2 contribute to the decay amplitudes (|4j) and (^) with different kinematic 
factors, we find that globally over all phase space (but with q 2 > 1 GeV 2 ) the kinematic 
factors are approximately of equal weight which makes it useful to compare the overall 
magnitudes of the form factors. We find ^> F^ 1 and F^ 3> For the b — > dq'q' 
amplitude we find that F± is dominant ~ 

The individual contributions \K^ d Kibfi(xi, q 2 )\, (i = u,c,t), to F^ are shown in Fig. ||. 
These magnitudes are the same for the b — > d + g transition. The c quark is the largest 
contributor. The weak phase from the CKM matrix is very small but this contribution 
carries a strong phase for g 2 > 4m 2 . This strong phase is the same for the b — > d + g 
transition. The u quark contribution has a weak CP-violating phase e~ l5vi ~ e~ i7 (Particle 
Data Group notation) and also a strong phase that is common to the b —>■ d + g transition. 
The t quark contribution is negligible. 

These individual amplitudes add to make F± and, because u and c make significant 
contributions, the phase of F± differs for the b — > d + g and b — > d + g transitions. The 
phase difference, which can be called the net CP-violating phase, is not negligible but has 
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no simple relationship with any particular angle of the unitary triangle. With S13 = ir/2 
and S13 = 0.0035, we show this phase in Fig. |3|. 

Because of the presence of both strong and weak phases the magnitudes of F± are also 
different for the b and b decays. Processes like b — > dss and b — > dss are expected to be 
penguin dominated and F± dominates all the other form factors. The decay rates dT/dq 2 
calculated from (gj) are shown in Fig. |j. The cc threshold cusp is clearly exhibited and 
CP violation is manifest. The difference of the decay rates can easily be shown to be 
proportional to the Jarlskog factor |25| ^s[K u bK* d K C( iK* b } = C12C13C23S12S13S23 sin 613. Since 



this factor basically controls the magnitude of the asymmetry, the modification with different 
choices of S13 and S13 (the least known elements of the CKM matrix) can be assessed. The 
asymmetry is large because the sum of the decay rates is also small. 

Turning to the b — *> sq'q' transition, we again find that F± 3> F^, F% 3> F 2 L and the F^ 
amplitude to be dominant. The individual contributions from u, c and t are shown in Fig. 
|5]. The c quark contribution in this case greatly outweighs that of the u and t quarks and 
since its contribution is so large and has almost zero weak phase, the weak phase on F^ is 
very small. Processes like b — > sdd and b — > sdd are expected to be penguin dominated and 
these lowest order calculations give the decay rates shown in Fig. |^. 

III. THE GLUINO PENGUIN IN THE MSSM 

In the MSSM there are contributions to from the two gluino exchange diagrams I and 
II (Fig. 0) corresponding to the gluon line attached respectively to the gluino and d squark 
lines. The MSSM penguin amplitudes have the form 

vn SM (g 2 ) = E Ajfo {C 2 (G)A l ABlM 

j 

+ [-C 2 (G) + 2C 2 (R)}A%^}P A (32) 

where (A, B) are chirality indices, C 2 (G) = 3 and C 2 (R) = E a T a T a = 4/3 are SU(3) 
Casimir invariants and j = 1, . . . , 6 labels the d squark mass eigenstates. The coefficient 
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\bq = 9 S rrjq * rrjb 
ABj - 4m 2 V dA V dB 



(33) 



describes the rotation from the down-diagonal interaction states to the d mass eigenstates 
at the d — d — g vertices. The matrices Vj L and Vg R are obtained from the (6 x 6) matrix 
= (Kit,? VdR) T which diagonalises the d mass 2 matrix 



Ml 



{Mj) LL {Ml) 



\ 



mu (M?) 



LR 



RR 



(34) 



Placing the external quarks on mass shell converts V^f SM (g 2 ) into the same general form 
(|13|) as for V (q 2 ) so that the MSSM form factors can also be obtained from (|14T) -( |16|) . 
For the LL MSSM penguin we find, after neglecting terms of order m 2 /m~ and m 2 /m~, 



a 



MSSM/ 2 



Zj(x,y) 



pMSSM {g2) = ^ ^ I dx 



-C 2 (^) + 2^)] 2 ^-f 1 -^ } ) 

^ (G > z^) 



and 



£ Al%- / dx / dy{C 2 (G) 



l-a; 



Zj(x,y) 



xy 



o Zj[x,y) 



where 



(35) 



(36) 



(37) 



(38) 



Zj(x,y) = 1 — x + Xj + g 2 [xy + y(y — l)]/m 



Z'Ax,y) = x + Xj(l -x) + q 2 [xy + y(y - l)]/m| 



(39) 
(40) 
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with Xj = ra| 

As g 2 <C m| we can set g 2 = in (|39| ) and ( ^0|) to get the LL penguin contributions 

F?(0) = £ A b L q Lj [C 2 (G)A(x,) + C*{R)B{x s )\, (41) 
j 

if (0) = 0, (42) 

— F 2 L (0) = —F 2 R 
m q m b 

M 



E Al%-[G 2 (G)C(%) - C 2 (itOL>(%)] (43) 



where 



A(x) = , 1 rj \3 -9x + 9x 2 - 3x 3 
K J 6(l-x) 4L 

+ (1 -3a; 2 + 2a; 3 )lnx], (44) 

B(x) = ~ ,. [11 - 18s + 9x 2 -2x 3 + 61na;], (45) 
18(1 — x) 4 

C(x) = 4 ^~ 1 ^ 3 [l-x 2 + 2xlna;], (46) 

D(x) = ~ , 4 [2 + 3x - 6x 2 + x 3 + 6s In a:]. (47) 

For the LR MSSM penguin 

x{C 2 (G) 



( 9 2 ) = -^E A L% / rfx/"^ 

m a ~" J Jo Jo 
x + 2y-l 



Zj(x,y) 

+[-C 2 (G) + 2C 2 (^] a; + 2y " 1 } (48) 

6j\x, y) 

(3 MSSM (q 2 ) = (49) 



and 



o 7o z j(x,y) 
+ [-C 2 (G) + 2C 2 (R)]^-^}, (50) 

if (g 2 ) = 0. (51) 

Again we can set q 2 = to obtain for the LR penguin contributions 
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Ff(0) = F?(0) = F*(0) = (52) 



and 



F 2 L (Q) =^E AS Rj [C 2 (G)E(x j ) - AC 2 (R)C(x 3 )} (53) 



with 



E(x) = — — ^— -[1 - x + xlnx}. (54) 
[I — x) 

The RR and RL penguins are obtained from the above by the replacements Allj ~~ > A^yy 
and A LRj — > A RLj respectively together with (-m g a MSSM ) <-> (m b( 5 MSSM ) and Fh^ <-> -F(f )2 )- 
The total g 2 = MSSM form factors are therefore 

^i L (0) = E Al%-[C 2 (G)A(%) + C 2 (it!)£(%)], (55) 
j 

F*(0) = £ A^.[C 2 (G)A(%) + 0,(^(5*)], (56) 

3 

^ 2 L (0)=E{K Al% +m 6 A^.] 

j 

x[a 2 (G)C(%)-C 2 ( J R) J D(%)] 

+m § Agy[C a (G)^(%) - 4C 2 (i?)C(£,)]}, (57) 

^(0)=E{K A l% +^ Agy] 

x[C 2 (G)C(^)-C 2 (i?) J D(x J )] 

+m § A^.[C 2 (G)£(%) - 4C 2 (i?)C(* i )]}. (58) 

The results for the F 2 L,i ^(0) MSSM form factors agree with those of |p6|| . However for 
the Fi (0) form factors, whereas our C 2 (i2) term is the same as that of |26] and [|TJ, 
the A(x) function occurring in the C 2 (G) term differs from that of by —(1 — a?)~ 2 lnx 
and bears little resemblance to the F(x) function of JTJ. Note though that our result for 
C 2 (G)A(x) + C 2 (R)B(x) is the same as the function Pp — |Pb given in f27jj . 

The MSSM calculations are described in ||. Two-loop MSSM RGEs were used for the 
gauge and Yukawa couplings and one-loop MSSM RGEs for the other SUSY parameters. 
Full flavour dependence was included in the running, with one- loop QCD and stop/gluino 

13 



corrections to the physical top mass from [28]. The unification scale boundary conditions 
were a universal scalar mass m , universal gaugino mass m 1( / 2 and a universal soft SUSY- 
breaking trilinear scalar coupling A. After minimisation at the scale m t of the full one- 
loop Higgs effective potential, which included all contributions from the matter and gauge 
sectors, we are left with a four-dimensional parameter space {m Q , 9711/2, A, tan/?}, where 
tan/? = V2/V1 is the ratio of the vacuum expectation values of the two Higgs fields, together 
with the sign of the coupling \i between the two Higgs fields. The physical Higgs masses 



were determined using the approximation to the RG- improved Higgs masses described in p9 
Mass eigenvalues and diagonalisation matrices for the d squarks were generated for a selection 
of data sets in the parameter space 150 < m < 1150; 150 < rrii/2 < 1150; 150 < \A\ < 1150 
(units of GeV) and 2 < tan/? < 48 which satisfied current experimental constraints (see 
and yielded a neutralino as the lightest supersymmetric particle. We also imposed 



the condition that the Standard Model like minimum be the global one as has become 



customary pi |. However it should be noted that, as pointed out in [32], this traditional 
condition is not sufficient to avoid cosmological problems. For this one should employ the 
slightly more restrictive condition in . 



The allowed values of A become more restricted by unphysical (charge and colour break- 
ing) minima as tan/3 increases from its fixed point value of tan/3 rs 1.5 |32 |. The avoidance 



of unphysical minima gives a bound of mo/wi/2 ^ 1 at the low fixed point which drops away 
to about 0.4 at intermediate values of tan/3. However the minimum bound on m /m 1 / 2 is 



for A ~ m and it increases quadratically in A away from this value ||32|| , so that effectively 
0.5 < A < 1.5mo at intermediate tan/3 values. Data sets for negative A were therefore more 
restricted in this region regardless of the sign of //, with all but those near mo producing 
colour breaking minima. Near the high tan (3 fixed point, where the bottom Yukawa coupling 



is large, the analysis of is 110 longer valid and the parameter space becomes once again 
less restricted here. Negative and quite large values of A are allowed (and even favoured) 
over positive ones in this region. 

Finally we should add that additional and probably very restrictive constraints on mo 
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especially at low tan (3 come from the need to avoid neutralino dark matter overdosing the 



universe. This was examined recently in |33|] but has not been included in our analysis here. 

The magnitudes of the MSSM form factors satisfy \F^\ > \F^\ > |F 2 L | ^> \F^\ for all 
regions of the allowed parameter space apart from the narrow region tan/5 = 2, mi/2 = 150 
and m ^ 1000 where \Fj | is slightly smaller than \F^\- Outside this region the ratio 
l-P^I/l-F^I exceeds unity and increases strongly with tan/3. For tan/3 = 2, the ratio ranges 
from Pd 2 for = 250 and niQ > 1000 to ~ 9 for low (mo, trixj-i) = (150,250). For 

higher tan (3, \F^\ becomes more dominant, the ratio increasing to 24-28 for tan/? = 10 and 
200-225 for tan f3 = 48. The relative sizes of the form factors are due to both the mixing 
coefficients A^ Bj - ( |33"D and the functions A, B, C, D and E of the variable Xj = m~ /m|. If 
the j dependence of Xj is neglected, the quantities A b ^ B = | J2j ^ABj I satisfy A b ^ L > A b RR > 
A b £ R > A b RL and this accounts in the main for the relative sizes of the form factors. The 
large values of the form factors at high tan (3 are due to an interplay of two factors; (i) the 
light gluino mass (rrig w 360 GeV) associated with = 150 and (ii) a gluino mass lying 
within the range of d masses such that the variable Xj is close to unity for several values of 
j- 

The result that F% is the largest MSSM form factor indicates that, in contrast to the 
SM, the magnetic dipole transition dominates the b decay process in the MSSM. To com- 
pare with the SM, we note that the ratio of the largest MSSM and SM form factors is 
|F 2 i? (MSSM)|/|FL(SM)(q 2 = 0)| < 0.4 GeV. 

The phases of the MSSM form factors change very little over the allowed parameter 
space. The phases of F± and F^ L ' R ' are independent of the sign of A and, for /i < 0, are 
approximately equal at ~ —2.8 for b — > d and w —0.016 for b — > s. For fi > the phases of 
F 2 ' are shifted by n. The phase of 

F R 

varies a little with mo and and depends on 
the sign of A, being approximately that of F^ for A > and shifted by 7r from that of F± 
for A < 0. These MSSM phases for // < are comparable to the corresponding SM phases 
so that the magnitude of the phase difference between the dominant MSSM form factor F 2 R 
and that of the SM form factor F^(q 2 = 0) is ~ 0.4 for b — > d and ~ 0.01 for b — > s for 
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fi < and ^2.7 and « 3.1 for /i > 0. Hence, after allowance for the negative sign in (p3|). 
we conclude that the superpenguins and ordinary penguins have the same sign for \i > 
and opposite sign for /i < 0. 



IV. SUSY EFFECTS IN b —* qq'q' 

One albeit crude measure of the effects of SUSY in the decays b — > gg'g' is the relative 
size of the integrated decay rates for the MSSM, taken in isolation, and for the SM. In 
computing these decay rates from @ we impose the low q 2 cutoff q 2 > 1 GeV to avoid 
non-perturbative long distance effects. 

The largest effects of SUSY on the decay rates occur for high tan (3 and low (mo, mi/ 2 ). 
For tan/? = 48 and A = -300 the ratio r Peng (MSSM)/r Peng (SM) has a maximum value at 
(m , m 1/2 ) = (275, 150) of w 0.10 (6 -> d) and « 0.085 (b -> s) for /i > and « 0.09 (6 -> d) 
and ~ 0.08 (6 — > s) for /x < 0. The ratio exceeds 10~ 2 at tan/3 = 48 for (m ,mi/ 2 ) 
ranging from (225, 150) to (275, 225). However, for lower tan/? this ratio has a much smaller 
maximum value; for tan/5 = 10 it is 1 x 10~ 4 at (mo,?ni/2) = (275,225) and 3 x 10~ 4 at 
(550, 150) for tan (3 = 2. The ratio decreases rapidly for large values of due mainly to the 
increase of the gluino mass in ( |33"1) from m 360 at m^j = 150 to ~ 1875 at m^j = 850. These 
findings for low and medium tan /3 differ from the earlier estimates U that the SUSY and 
SM contributions to T(b — > sq'q') were of comparable size. However, these early estimates 
were based on the assumption that b — > sq'q' could be described solely by the LL penguin 
form factors Fi ' R \o) and our study shows F 2 to be the dominant form factor. 

In SUSY there is also a contribution to b — > qq'q' from the box diagrams of Fig. ^| for 
which the amplitude is |],|23 



• 2 

mB ° X = ^2 [u q (p q )T a l,P L u b (p b )] {J l LL [u ql {jp ql )^T a P L v ql {j> ql )] 

+Jl R [{%'(P q >h'*T a P R V q -,{pq,))} + [u q (p q )^P L U b ( Pb )} 

x{4l [uq'M^PLVfM} + Jl R [u q >(p q >)YP R v q >(p q >)}} + {L^R) (59) 



16 



where 



Jab = E E U^Vi) Vj& * vg (60) 



with 



7 2 

"A = g ft) - g 2/i), (61) 

1 7 

^2 = -g g(Xj,Vi) + g /^J, j/t), (62) 

2 4 

^3 = g 2/i) + g /(%, Jfi). (63) 

Here & = m~./m~, the squark q is w for q' — u and d for q' = d and the box functions are 



11 x 

g(x,y) = [ - + ( -) 2 lnx- (x ^ y)] (64) 

x — y y — 1 x — 1 

11 x 

f(x,y) = [ - + -. —r\nx-(x^y)\. (65) 

y — x y — 1 (x — I)" 1 

For the allowed regions of parameter space the box amplitudes satisfy \Jl R \ > 3(| J| L |, 
\ j lr\) » (\JrlI \ j rr\i \ j rl\i \ j rr\) a P art from in the region tan/5 = 2, m 1/2 = 150 
and mo > 650 where | J^lI becomes slightly larger than The four largest box ampli- 

tudes are generally of the same order as the MSSM penguin amplitudes (Ff, F 2 {L ' R) ); the 
remaining four are negligible, being smaller by a factor of at least 10 5 and comparable to 
the MSSM penguin Ff. For the regions where r MSSM /T SM > 10 -3 , the ratio of the largest 
MSSM box and penguin amplitudes \Jlr\/\F 2 r \ * s small, varying from only 8 x 10~ 3 for the 
parameters tan/3 = 48, (mcmift) = (275, 150) which produce the maximum SUSY effects 
to a maximum of 0.01. The ratio does increase to ~ 0.4 for low tan (3 (=2), low mi/2 (= 150) 
and high m > 1000 but in these regions the SUSY effects are negligible. Hence, for the 
MSSM data sets for which the SUSY penguin effects are largest, the SUSY box amplitudes 
can be neglected in calculating the decay rates. 

The differential decay rates dT Pens /dq 2 (|6[) for the SM and for the combined effects of 
the SM and the MSSM for the MSSM data set A = -300, \i > 0, tan/3 = 48, (m , m 1/2 ) = 
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(275, 150) which maximises the SUSY effects show (see Fig. f|and the SUSY enhancement 
of the decay rates to be significant for most of the range of q 2 values. The partial rate CP 
asymmetries Acp(q 2 ), defined in ([TOD , reveal the presence of SUSY for q 2 < Am 2 . For these 
values of q 2 the SM CP asymmetries of ~ 25% for b — > d and « 1.5% for b —>■ s are reduced 
to ~ 20% and ~ 1.2% respectively when the MSSM contributions are included. 



V. DISCUSSION AND CONCLUSIONS 

We have calculated, from first principles, both the SM and MSSM penguins that con- 
tribute to the decays of the b and b quarks. For the MSSM in particular there are discrepan- 
cies to be found in the literature p6|fl as to the correct formulae for the Fi' R ' form factors. 



Our results for these form factors differ slightly from those of [26] but our other results agree 



with p6 



Because of the presence of strong phases in the contributions to the SM penguin ampli- 
tudes from u and c quarks we find that the decay rates for b — > dq'q' is significantly different 
to the rate b —>■ dq'q' even for quarks in isolation (see Fig. [(]). 

The SUSY enhancement of the gluon-mediated exclusive hadronic b decays within the 
constrained MSSM model can be at the several percent level in certain regions of the 
(A, tan/3, m , mx/2) parameter space. In these regions the SUSY penguin processes 
dominate the SUSY box processes with the consequence that the b decays in the MSSM 
are driven by the magnetic dipole transition rather than the electric dipole transition of the 
SM. 

QCD corrections arising from renormalisation of the present short distance results down 
from the electroweak scale to the scale rrib are not likely to alter the finding that the magnetic 
amplitude is dominant in the MSSM as the QCD induced mixing effects [|T2"| , |3lf produce an 



enhancement of the magnetic dipole operators in the AB = 1 effective Hamiltonian relative 
to the current-current penguin operators associated with the electric dipole amplitude. Fur- 
thermore, Gerard and Hou 0] have noted that the result fl2"P| ) for the SM form factor F^(q 2 ) 
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already contains the dominant part of the QCD corrections for the current-current penguin 
operators and, therefore, that the main effects of QCD corrections will be the renormalisa- 
tion of the strong coupling constant from a s (Mw) to a s (mi,). This would have the effect of 
increasing the penguin decay rates of the SM by the factor rj = a s (m b ) /a s (M w ) ps 1.84 and 
also increasing the MSSM penguin amplitudes relative to those of the SM. 

Detection of new physics in the hadronic decay amplitudes of the b quark through a study 
of deviations from the predictions of the SM in the patterns of CP violation in B d decays 
is complicated, on the one hand by the interplay between the cumulative effects in the SM 
of the g 2 -dependent strong phases in and the weak CKM phases from the contributing 
u, c and t quarks and, on the other hand, by MSSM phases comparable in magnitude to 
the SM weak phases and which can give constructive or destructive interference depending 
upon the details of the soft SUSY-breaking mechanism. 
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FIGURES 




FIG. 1. SM gluon penguins with W and scalar exchange. 
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FIG. 2. Contributions to SM Fi(q 2 ) forb^d + g from u (solid line), c (dashed line) and t 
(dotted line) quarks. 
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FIG. 3. Phase of SM F^{q 2 ) for b -> d + g (solid line) and b -> d + 5 (dashed line) for a CKM 
phase of tt/2. The dotted line shows the CP-violating phase difference. 
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FIG. 4. Differential SM decay rates for b — > dss (solid line) and b — > Jss (dashed line). The 
corresponding results for the combined effects of the SM and MSSM are given by the dotted and 
dot-dashed lines respectively. 
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FIG. 5. Contributions to SM Fi(q 2 ) for b — ► s + g from u (solid line), c (dashed line) and £ 
(dotted line) quarks. The contribution from the c quark has been scaled down by a factor of 10. 
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FIG. 6. Differential SM decay rates for b — > sdd (solid line) and 6 — > sd<i (dashed line). The 
corresponding results for the combined effects of the SM and MSSM are given by the dotted and 
dot-dashed lines respectively. 
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FIG. 7. MSSM gluino penguins. 
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FIG. 8. MSSM box contributions to b -> qq'q'. 
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